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Abstract 

The nonlinear ballistic conductance of three-dimensional quantum microcon- 
strictions, which contain magnetic impurities, is investigated. The nonlinear 
part of the conductance, which is due to the interaction of electrons with 
magnetic impurities is obtained. The analytical results have been analyzed 
numerically. It is shown that intensity of the Kondo anomaly in the conduc- 
tance as the function of the applied voltage depends on the diameter of the 
constriction and the positions of impurities. 

The impurity-electron interaction in Kondo systems can be effectively studied by using 
the point contacts (PC). In first measurements of the differential PC resistance R{V) in 
metals with magnetic impurities the zero-bias Kondo anomaly had been observed 
These experiments were explained by quasiclassical theory of Kondo effect in PC's [Q. It 
was shown that in second-order Born approximation the magnetic impurity contribution to 
the PC resistance includes the logarithmic dependence R{V) ~ In {V) for eV ^ Tk and 
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the saturation for eV <^ Tk (Tk is the Kondo temperature, V is the voltage apphed to the 
PC). In accordance to the theory [^], the nonhnear correction to the balhstic PC resistance 
is proportional to the contact diameter. But in the experiments the size dependence 
of the PC current was not investigated due to the limited range of contact diameters, which 
were accessible. 

The development of the technique of mechanically controllable break junctions (MCBJ) 
has made it possible to create the stable PC's, with the diameter adjustable over broad 
range, down to a single atom |^J^. In the MCBJ experiments authors had studied the 
resistance of ultrasmall contacts with magnetic impurities as function of the PC diameter d. 
In the contrast to the prediction of the quasiclassical theory |^ Yanson at el. |^,|| observed 
that Kondo scattering contribution to the contact resistance is nearly independent from the 
contact diameter d for small d. Such behavior authors had explained by the increasing 
of Kondo impurity scattering cross-section with decreasing of contact size. 

In theoretical works |p it was shown that in very small contacts the discreteness of 
impurity positions must be taken into account, and experiments may be explained by 
the "classical" mesoscopic effect of the dependence of the point contact conductance on 
the spatial distribution of impurities. This effect is essential in the "short" contacts and 
in the quasiclassical approximation it disappears with the increasing of the contact length. 
Zarand and Udvardi had considered the contact in the in the form of a long channel and 
suggested that the Kondo temperature is changed due to the strong the local density 
of states fluctuations generated by the reflections of conduction electrons at the surface of 
the contact. As a result of that, the effective cross section of electrons has the maximum, 
if the position of the impurity inside the contact corresponds to the maximum in the local 
electron density of states. But the mesoscopic effect of the spatial distribution of impurities 
in quantum contacts was not analyzed in the paper . 



In ultrasmall contacts the quantum phenomena, which are known as quantum size effect, 
occur. The effect of the 2e^ /h conductance quantization has been observed in experiments on 
contacts in the two-dimensional electron gas [|ll] , |T2| and in the ultrasmall three-dimensional 



constrictions, which is created by using the scanning tunnel microscopy |1T3|,[T^ and mechan 



ically controllable break junctions [|I^. The defects produce the backscattering of electrons, 
and thus break the quantization of the conductance. From the other hand, the impurities 
situated inside the quantum microconstriction produce the nonlinear dependence of the con- 



ductance on the applied voltage |TB|. This dependence is the result of the interference of 
electron waves reflected by these defects [^,0]. 

In this paper we present the theoretical solution of this problem for the conductance 
of a quantum microconstriction in the form of the long ballistic channel , which contains 
single magnetic impurities. The study is made of the first and second order corrections to 
the conductance of the ballistic microconstriction in the Born approximation. The effect of 
impurity positions is taken into account. Within the framework of the model of the long 
channel the quantum formula for the conductance G is obtained. By using the model of the 
cylindrical microconstriction, the nonlinear conductance as a function of voltage V and the 
width of constriction d is analyzed numerically for different positions of a single impurity. 

Let us consider the quantum microconstriction in the form of a long and perfectly clean 
channel with smooth boundaries and a diameter d comparable with the Fermi wave length 



Xp = h/y^Q/mep, where ep is the Fermi energy. We assume that this channel is smoothly 
(over Fermi length scale) connected with a bulk metal banks. As it was shown [^,^], in 
such constriction in the zeroth approximation on the adiabatic parameter | V(i| ^ 1 accurate 
quantization can be obtained. The corrections to the tunneling and reflection coefficients 
of electrons due to deviation from the adiabatic constriction are exponentially small, except 
near the points where the modes are switched on and off . 



When a voltage V is applied to the constriction, a net current start to flow. In the limit 
V —>■ 0, the ballistic conductance of the quantum microconstriction is given by the formula 

where fp is the Fermi function, ep is the minimal energy of the transverse electron mode, f3 is 
the full set of transverse discrete quantum numbers. The ballistic quantum PC displays the 



specific nonlinear properties, such as the conductance jumps jh. For the two dimensional 



PC these effects was considered in the papers [ p3| , p^ . The aim of this study is to analyze 
the zero bias Kondo minimum in the PC conductance. We assume that the bias eV is 
much smaller not only the Fermi energy ep-, but also the distances between the energies 
Bp of quantum modes. In this case the effect of the influence of the applied bias to the 
transmission is negligibly small. 

Impurities and defects scatter the electrons that leads to the decreasing of the transmis- 
sion probability. In accordance with the standard procedure p5| , p6| the decreasing of the 
electrical current A/ due to the electron-impurity interaction connects with the velocity of 
the energy E dissipation by the relation: 



dt 



d{H,) 
dt ' 



The Hamiltonian of the electrons H contains the following terms:. 



(2) 



H = Ho + Hi + H, 



inti 



(3) 



where 



is Hamiltonian of free electrons, 



Ho =^ ekcl^Cko 



T 



(4) 



(5) 



describes the interaction of electrons with electric field. The Hamiltonian of the interaction 
of electrons with magnetic impurities Hint can be written as 



Hint = J2 ^JAk' [Sz (cl^Ckr - cl^Cki) + S cl^Cki + +S+cl,^Ck^ 



(6) 



j,k,k' 
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Here the operator c^^ (cka) creates (annihilates) a conduction electron with spin a, wave 
function (pk, and energy Sk] S denotes the spin of impurity; is the electron velocity along 
the channel; Jj ^ j^.' is the matrix element of the exchange interaction of electron with impurity 
in the point r^ ; ka is the full set of quantum numbers; 

Jj,k,k' ^ J Ji^^ ^jWk{r)^l'{r). (7) 

The electron wave functions and eigenvalues in the long channel in the adiabatic approxi- 
mation are 



<^fe(r) = ^/3(R) exp (^^PzZ^ ; (8) 

^' = ^'' + ^' 

where k — {l3,Pz) , P is the set of discrete transverse quantum numbers; Pz is the momentum 
of an electron along the contact axis; r/ig is an electron mass; r = (R,,z) , R is a coordinate 
in the plain, perpendicular to the z axis. 

Differentiating (Hi) over the time t we obtain the equation for the changing A7 of the 
current as a result of the interaction of electrons with magnetic impurities: 



VAI = ^{[H,{t),Hint{t)]), (10) 
in 

where 

{...)= Trip it)...). (11) 

All operators are in the representation of interaction. 

The statistical operator p (t) satisfies to equation 



lh-^ = [Hint {t),p{t)], 
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(12) 



which can be solved using the perturbation theory: 

t t t' 



(13) 



-oo — oo 



By means of Eq.(13) the changing in the electric current due to magnetic impurities can 
be determined 



A/ = /i + /2 + ... = (14) 
^ J dt'Tr {pom,H,„,{t)],Hi^t{t')])- 



—oo 

t' t 

^ • dt" / dt'Tr (po [[[i^i, i^^n^(^)] , H^^tit')] , HUt")]) + - 



-oo — oo 



After the simple, but cumbersome calculations we find the first and second order correc- 
tions to the PC current 



h = ~^^(^ + ^) XI X] ^^"^^^ ^^"^ ~ '"^n)ifm - fn) S {Sn - 8^) J j,n,m ^i,m,n] 

n.m i.j 

(15) 



I2 = ys(s + 1) X X^*^^" - ^^gn v,J (16) 

n,m,k 

[5(£„-£fc)Pr h5(£m-£jfc)Pr 



[2/n(/fe — /m) + (/m " fk)], 

where = fp {sn + ^signvz) ■ The first addition 7i to the PC current describes the small 
spin-depended correction (of the order (J/ep)^ ) to the changing of the current due to 
the usual scattering. The second addition I2 is also small too, but contains the Kondo 



logarithmic dependence on the voltage, and it is most important for the analysis of the 
nonlinear conductance of constrictions with magnetic impurities. 

The expressions (15) and (16) can be further simplified in the case of 5— potential of 
impurities 

J(r) = J5(r) (17) 
In this case the addition I2 to the ballistic current has the form: 



2Jhie_ 

h 



n,m,k 



(18) 



[5{£n - £k) Pr ^ V5{em- Sk) Pr — - 



£k 



[2/n(/A; — fm) + {fm — fk)], 

As it follows from the Eqs.(16), (18), the current I2 depends from the positions of im- 
purities. Two effects influence by value I2 : the effect of quantum interference of scattered 
electron waves, which depends from the distances between impurities, and effect of the elec- 
tron density of states in the points, where the impurities are situated. The nonlinear part 
of the conductance can be easy obtained after differentiation the Eq.l8 over the voltage 
G2 = dl2/dV. In the case of a single impurity and at zero temperature T = this equa- 
tion can be analytically integrated over momentum and the conductance G2 takes the 
following form: 



ne ml 



a, 13,"/ >t=± 



Pa Pp P-y 



(19) 



In 



Pi Pi [Pa 



p\ +p\ 



+ (1 - 5o,p) In 



Pa Pp - Pa Pp 



Pa Pp + Pa Pp 



+ 



^ap In 



P 



(^) 



p. 



where 

p(f ) = ^2m, (^eF±^-eo^, (20) 

and the transverse parts of the wavefunction ipai^) and the electron energy Sa are defined 
by Eqs. (8), (9). 

Carrying out the numerical calculations we use the free electron model of a point contact 
consisting of two infinite half-spaces connected by a long ballistic cylinder of a radius R and 
a length L (Fig.l). In a limit L — oo the electron wave functions (fk (r) and eigenstates Sk 
can be written as 

(r) = /-^ ^ , -Jm (imn-^) exp (im^ + ^Pzz) ] (21) 

V^Jm+lilmn) ^ \ Tl J 

and cylindrical coordinates r = (p, ip, z) with the axis z along the channel axis have been 
used. Here k = {n, m,pz) are the quantum numbers, Q = ttR'^L is the volume of the channel, 
7mn are the n-th zero of the Bessel function J^. Because the degeneration of the electron 
energy on azimuthal quantum number m ( as a result of the symmetry of the model), 
quantum modes with ±m give the same contribution to the conductance. In this model the 
ballistic conductance (1) has not only steps Go, but also steps 2Gq . 



In Fig. 2 the dependence of the nonlinear conductance on the applied bias is shown 
for the different positions of a single magnetic impurity inside the channel. The results 
obtained confirm that the nonlinear effect is strongly depend on the position of impurity. 
If the impurity is situated near the surface of the constriction r = R, where the square 
module of the electron wave function is small, its influence to the conductivity is negligible. 
This conclusion is conflrmed by the calculations of the dependence G2 on the position of 
the impurity for different number of quantum modes (Fig. 3). Results indicate that the 
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mesoscopic effect of the impurity position is more essential for ultrasmall contacts, which 
contain only few conducting modes, and has a maximum. The similar results is obtained 
for the dependence of G2 on the radius R of the constriction (Figs. 4, 5). In the single-mode 
constriction (Fig.4) the conductance G2 displays much more stronger dependence on i?, than 
in the contact with five conducting modes (Fig.5). 

Thus, we have shown that in the long quantum microconstrictions the spatial distribution 
of magnetic impurities influences to the nonlinear dependence of the conductance on the 
applied voltage. This mesoscopic effect is due to the strong dependence the amplitude of 
an electron scattering on the positions of impurities. As a result of the reflection from the 
boundaries of the constriction the electron wave functions, which correspond to the flnite 
electron motion in the transverse to the contact axis direction, are the standing waves. If 
the impurity is situated near the point, in which the electron wave function is equal to zero 
(near the surface of the constriction or, for quantum modes with numbers n > 1, in some 
points inside), its scattering of electrons is small. The fact the amphtude of the Kondo 
minimum of the conductance of the quantum contact display the mesoscopic effect of the 
dependence on the positions of single impurities. This effect is most important in the case, 
when only few quantum modes are responsible on the conductivity of the constriction. 
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Figure captions. 

Fig. 1. Schematic representation of a ballistic microconstriction in the form of a long 
channel, adiabatically connected to large metallic reservoirs. Magnetic impurities inside the 
constriction are shown. 

Fig. 2. The voltage dependence of the nonlinear part of conductance G2 (19) from the 
distance of the impurity from the contact axis {2'nR = 5.2Xf; T = 0; 1 - 27rp = l.bXp ; 2 
-27rp = 2.5Xf ; 3 -27rp = 3.0Af ; 4 - 27rp = 3.5Af ). 

Fig. 3. The dependence of G2 (19) on the position of the impurity for the different 
quantum modes in the constriction {V = 0.02sf', T = 0; 1 - one mode {2ttR = SXp)', 2 - 
three modes(27rit! = 4Ai?); 3 - five modes {27rR — 5.3Ai;'); 4 - six modes {27rR — QXp))- 

Fig. 4. The dependence of G2 (19) on the radius of the constriction for the single mode 
channel and different positions of the impurity {V = 0.02£f;T = 0;1 - 27rp = O.SA^;' ; 2 
-27rp = I.OXf ; 3 -27rp = 1.5Af ; 4 - 27rp = 2.0Af ) 

Fig. 5. The dependence of G2 on the radius for the microconstriction with five quantum 
modes and different positions of the impurity {V = 0.02sf;T = 0;1 - 27r/9 = O.bXp ; 2 
-27rp = 1.5Xf ; 3 -27rp = 2.5Xf ; 4 - 27rp = 4.5Af ) 
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